Stochastic Model for Tumor Control Probability: Effects of Cell Cycle and (A)symmetric Proliferation

Estimating the required dose in radiotherapy is of crucial importance since the administrated dose should be sufficient to eradicate the tumor and at the same time should inflict minimal damage on normal cells. The probability that a given dose and schedule of ionizing radiation eradicates all the tumor cells in a given tissue is called the tumor control probability (TCP), and is often used to compare various treatment strategies used in radiation therapy. In this paper, we aim to investigate the effects of including cell-cycle phase on the TCP by analyzing a stochastic model of a tumor comprised of actively dividing cells and quiescent cells with different radiation sensitivities. We derive an exact phase-diagram for the steady-state TCP of the model and show that at high, clinically-relevant doses of radiation, the distinction between active and quiescent tumor cells (i.e. accounting for cell-cycle effects) becomes of negligible importance in terms of its effect on the TCP curve. However, for very low doses of radiation, these proportions become significant determinants of the TCP. Moreover, we use a novel numerical approach based on the method of characteristics for partial differential equations, validated by the Gillespie algorithm, to compute the TCP as a function of time. We observe that our results differ from the results in the literature using similar existing models, even though similar parameters values are used, and the reasons for this are discussed.Comment: 12 pages, 5 figure

Phenotypic heterogeneity in modeling cancer evolution

The unwelcome evolution of malignancy during cancer progression emerges through a selection process in a complex heterogeneous population structure. In the present work, we investigate evolutionary dynamics in a phenotypically heterogeneous population of stem cells (SCs) and their associated progenitors. The fate of a malignant mutation is determined not only by overall stem cell and differentiated cell growth rates but also differentiation and dedifferentiation rates. We investigate the effect of such a complex population structure on the evolution of malignant mutations. We derive exact analytic results for the fixation probability of a mutant arising in each of the subpopulations. The analytic results are in almost perfect agreement with the numerical simulations. Moreover, a condition for evolutionary advantage of a mutant cell versus the wild type population is given in the present study. We also show that microenvironment-induced plasticity in invading mutants leads to more aggressive mutants with higher fixation probability. Our model predicts that decreasing polarity between stem and differentiated cells turnover would raise the survivability of non-plastic mutants; while it would suppress the development of malignancy for plastic mutants. We discuss our model in the context of colorectal/intestinal cancer (at the epithelium). This novel mathematical framework can be applied more generally to a variety of problems concerning selection in heterogeneous populations, in other contexts such as population genetics, and ecology.Comment: 28 pages, 7 figures, 2 table

Modeling the Spatial Distribution of Chronic Tumor Hypoxia: Implications for Experimental and Clinical Studies

Tumor oxygenation status is considered one of the important prognostic markers in cancer since it strongly influences the response of cancer cells to various treatments; in particular, to radiation therapy. Thus, a proper and accurate assessment of tumor oxygen distribution before the treatment may highly affect the outcome of the treatment. The heterogeneous nature of tumor hypoxia, mainly influenced by the complex tumor microenvironment, often makes its quantification very difficult. The usual methods used to measure tumor hypoxia are biomarkers and the polarographic needle electrode. Although these techniques may provide an acceptable assessment of hypoxia, they are invasive and may not always give a spatial distribution of hypoxia, which is very useful for treatment planning. An alternative method to quantify the tumor hypoxia is to use theoretical simulations with the knowledge of tumor vasculature. The purpose of this paper is to model tumor hypoxia using a known spatial distribution of tumor vasculature obtained from image data, to analyze the accuracy of polarographic needle electrode measurements in quantifying hypoxia, to quantify the optimum number of measurements required to satisfactorily evaluate the tumor oxygenation status, and to study the effects of hypoxia on radiation response. Our results indicate that the model successfully generated an accurate oxygenation map for tumor cross-sections with known vascular distribution. The method developed here provides a way to estimate tumor hypoxia and provides guidance in planning accurate and effective therapeutic strategies and invasive estimation techniques. Our results agree with the previous findings that the needle electrode technique gives a good estimate of tumor hypoxia if the sampling is done in a uniform way with 5-6 tracks of 20–30 measurements each. Moreover, the analysis indicates that the accurate measurement of oxygen profile can be very useful in determining right radiation doses to the patients

Mathematical modelling and biomechanics of the brain

This monograph aims to provide a rigorous yet accessible presentation of some fundamental concepts used in modeling brain mechanics and give a glimpse of the insights and advances that have arisen as a result of the nascent interaction of the mathematical and neurosurgical sciences. It begins with some historical perspective and a brief synopsis of the biomedical/biological manifestations of the clinical conditions/diseases considered. Each chapter proceeds with a discussion of the various mathematical models of the problems considered, starting with the simplest models and proceeding to more complex models where necessary. A detailed list of relevant references is provided at the end of each chapter. With the beginning research student in mind, the chapters have been crafted to be as self-contained as possible while addressing different clinical conditions and diseases. The book is intended as a brief introduction to both theoreticians and experimentalists interested in brain mechanics, with directions and guidance for further reading, for those who wish to pursue particular topics in greater depth. It can also be used as a complementary textbook in a graduate level course for neuroscientists and neuroengineers

Modeling the Spatial Distribution of Chronic Tumor Hypoxia: Implications for Experimental and Clinical Studies

Tumor oxygenation status is considered one of the important prognostic markers in cancer since it strongly influences the response of cancer cells to various treatments; in particular, to radiation therapy. Thus, a proper and accurate assessment of tumor oxygen distribution before the treatment may highly affect the outcome of the treatment. The heterogeneous nature of tumor hypoxia, mainly influenced by the complex tumor microenvironment, often makes its quantification very difficult. The usual methods used to measure tumor hypoxia are biomarkers and the polarographic needle electrode. Although these techniques may provide an acceptable assessment of hypoxia, they are invasive and may not always give a spatial distribution of hypoxia, which is very useful for treatment planning. An alternative method to quantify the tumor hypoxia is to use theoretical simulations with the knowledge of tumor vasculature. The purpose of this paper is to model tumor hypoxia using a known spatial distribution of tumor vasculature obtained from image data, to analyze the accuracy of polarographic needle electrode measurements in quantifying hypoxia, to quantify the optimum number of measurements required to satisfactorily evaluate the tumor oxygenation status, and to study the effects of hypoxia on radiation response. Our results indicate that the model successfully generated an accurate oxygenation map for tumor cross-sections with known vascular distribution. The method developed here provides a way to estimate tumor hypoxia and provides guidance in planning accurate and effective therapeutic strategies and invasive estimation techniques. Our results agree with the previous findings that the needle electrode technique gives a good estimate of tumor hypoxia if the sampling is done in a uniform way with 5-6 tracks of 20–30 measurements each. Moreover, the analysis indicates that the accurate measurement of oxygen profile can be very useful in determining right radiation doses to the patients.Peer Reviewe